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Abstract 

We start by a review of the chronology of mathematical results on the Dirichlet- 
to-Neumann map which paved the way towards the physics of transformational 
acoustics. We then rederive the expression for the (anisotropic) density and bulk 
modulus appearing in the pressure wave equation written in the transformed co- 
ordinates. A spherical acoustic cloak consisting of an alternation of homogeneous 
isotropic concentric layers is further proposed based on the effective medium the- 
ory. This cloak is characterised by a low reflection and good efficiency over a large 
bandwidth for both near and far fields, which approximates the ideal cloak with a 
inhomogeneous and anisotropic distribution of material parameters. The latter suf- 
fers from singular material parameters on its inner surface. This singularity depends 
upon the sharpness of corners, if the cloak has an irregular boundary, e.g. a polyhe- 
dron cloak becomes more and more singular when the number of vertices increases if 
it is star shaped. We thus analyse the acoustic response of a non-singular spherical 
cloak designed by blowing up a small ball instead of a point, as proposed in [Kohn, 
Shen, Vogelius, Weinstein, Inverse Problems 24, 015016, 2008]. The multilayered 
approximation of this cloak requires less extreme densities (especially for the lowest 
bound). Finally, we investigate another type of non-singular cloaks, known as invis- 
ibility carpets [Li and Pcndry, Phys. Rev. Lett. 101, 203901, 2008], which mimic 
the reflection by a flat ground. 

Keywords: Mathematical methods in physics; Metamaterials; Anisotropic optical materials. Waves, 
Acousto-optical devices; Acousto-optical materials, Scattering, invisibility; Invisibility cloaks. 
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1 Introduction: Dirichlet-to-Neumann map versus 
cloaking 



Several mathematical studies of invisibility devices have emerged in the early eighties from a consideration 
of electric impedance tomography (EIT) models [1] [51 [31 13] . The aim of this technique is to probe the 
spatial variations in conductivity cr within some three-dimensional bounded region D, by applying a known 
static voltage u to the surface dfl bounding the region and recording the resulting current aVu. The 
equation satisfied by these quantities inside f2 is 



where n is the unit outward normal to the surface dft. In order for EIT to work, it must be possible to 
determine a from a knowledge of A^. If this can be done, then cloaking (in this situation) is impossible, 
even in theory. The question of whether or not the mapping can be used to determine the form of 
the conductivity is known as the Calderon problem 5j. It turns out [5] that the Dirichlet-to-Neumann 
map does determine a , but only under certain conditions: namely, that a must be known to be scalar- 
valued, positive and finite However, if some of these conditions are not met (for instance a is matrix 
valued), cloaking may be possible. The singularity and anisotropy of the material parameters required 
for such a conductivity cloak allow for the possibility that the Dirichlet-to-Neumann map does not 
uniquely determine the conductivity [8j. In 2006, Leonhardt proposed a design for a two-dimensional 
electromagnetic cloak which works in the geometric optics limit. Here, the propagation of light is governed 
by the Helmholtz equation, for which a similar uniqueness theorem has been established using tools of 
quasi-conformal maps [S]- Electromagnetic sensing in the quasistatic limit, which is closely related to 
cloaking, was reported by the mathematician Lassas back in 1997 The analysis of Lassas, based on 
variational principles (see also [10] and [IT] for the all-frequency case) uses the fact that mathematically, 
the trace map in ^ associated with the Helmholtz equation V • (crVu) — oj'^uis well defined and invertible 
from the Hilbert space F-i/2(f]) to the Hilbert space H^/^{dVt) provided that the conductivity cr is a 
uniformly elliptic symmetric-matrix- valued function and the spectral parameter oj^ avoids a discrete set 
of eigenvalues of the Neumann Laplacian (which can be interpreted as eigenfrequencies associated with 
almost-trapped eigenstates u |12j in the context of quantum mechanics): This comes from the fact that 
the aforementioned Helmholtz equation subject to the boundary condition du/dn — f forms a well-posed 
boundary value problem and < C^\\f\\H-^/^{dn} where C is a strictly positive constant depending 

upon fl (assumed to be smooth enough, which excludes the case of bad domains, such as star-shaped 
as numerically studied in section 6 of this paper). Actually, a singular change of variable resulting in 
perfect cloaking in the context of electric impedance tomography was introduced by Greenleaf et al in 
space dimension n > 3. This was a variation of an abstract two-dimensional example of Lassas, et al 
[14]. Kohn et al show in [^ that the situation is not significantly different when n = 2: perfect cloaking 
is also possible in space dimension two, using tools for the removability of singularities for harmonic 
functions. In the resonance domain, whereby the wavelength of the wave incident on the cloak is about 
the same size as its diameter, it is actually necessary to use vector Maxwell's equations to account for 
effects of polarization and this was done back in 2006 by Pendry, Schurig and Smith [T5] . Although this 
last paper did not contain a mathematically clean derivation of the problem, a subsequent experimental 
paper involving the same authors ^16j was convincingly demonstrating that cloaking is possible, at least 
for microwaves, and at a specific frequency, 8.5 GHz. A full mathematical account of the invisibility 
problem in electromagnetism can be found in a paper by Greenleaf, Kurylev, Lassas and Uhlmann |10j . 
which is written with a strong mathematical flavor. 

^It is generally assumed in the physics literature that cr is a smooth function, with regularity depending 
on the dimension of the problem. In dimension 2, however, mathematicians have established that no 
smoothness is required: only boundedness and positivity are needed [7]. 



V • (crVii) . 



(1) 



The current-voltage relationship provides a Dirichlet-to-Neumann map A^. on dfl: 




(2) 
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In [T], it is shown (with a credit therein to Luc Tartar) that the essence of indistinguishable media|^is 
encompassed in boundary measurements in a very simple way. We find it is worthwhile reproducing the 
main argument here: Let F denote the map from the original set of coordinates to the new one. If F{x) = 
X at dil, then the boundary measurements associated with a and d = DF{x)a{x)DF[x)'^ /det{DF{x)) := 
JaJ'^ /det{3) are identical, keeping in mind the coefficients Jj — djji/dxj of the Jacobian matrix are 
evaluated at point x ~ F~^{y). In other words, we are assured that ^a{f) — Ag-(/) for all / Dirichlet 
data of ([I]). 

The proof of this statement ^ is straightforward since if F{x) = a; at 951, then the change of variables 
F does not affect the Dirichlet data. If it{x) is bounded positive and definite, then the solution of ([T]) 
with Dirichlet data / is given by the minimum of the energy functional J^^ a{x)\I ■ Vxudx, so for any 
/, we have 

/ fk„fdx^ min / a{x)V -V xudx min / cr{y)V yU ■ V yudy = / fA^fdy. (3) 
Jan "=/ at an Jn «=/ at an Jn Jan 

Importantly, this means the Dirichelt-to-Neumann map can determine a at best up to a change a 
variables, as first noted by Kohn and Vogelius in 1984 Ij with a reference therein to Luc Tartar: and 
Act determine identical quadratic forms on Dirichlet data. This is the essence of cloaking. If, however, 
a is scalar- valued, positive, and finite, then the situation is very different: under some conditions on 
the regularity of fl, knowledge of the Dirichlet-to- Neumann map A^- determines an internal isotropic 
conductivity a uniquely, and cloaking breaks down. The salient consequence is that any cloaking device 
requires anisotropy, a fact known since the eighties, due to Sylvester and Uhlmann [H] for n > 3 in 1987 
and Nachman [19 for n = 2 in the 1995, but rediscovered and popularized twenty years afterwards in 
the context of invisibility cloaks by Pendry and his colleagues. 

In the present paper, we are actually interested in pressure waves propagating in fluids. Such a 
problem is of a different physical nature than electromagnetic waves interacting with invisibility cloaks 
(e.g. there are no polarization effects), but the associated governing equation (the scalar Helmholtz 
equation) is a simplified version of Maxwell's equations [ 10] . This might explain why acoustic cloaking 
has been touched upon by various photonics groups including |201 121j in two-dimensional settings and 
[231 US] in three-dimensional settings. Using tensors notations with Einstein's summation convention, 
one can write down the form-preserving transformations for the density and bulk modulus as: 

= detiJl) ' detiJl) ' ^ > 

with a indicating that it is a rank-2 tensor. This very useful formula is rederived using simple arguments 
in section 2. 

We emphasize that the formula Q associated with transformation based solutions to the conductivity 
and acoustic equations is a version of the form preserving transformation for the tensors of permittivity 
and permeability e^^ = {J^e'^'' ) / det{Jl) and /i'-' — {Jl^i'''' Jf) / det{Jl) appearing in Maxwell's equations 
in curvilinear coordinate systems, which were used by Pendry, Schurig, and Smith [TS] in order to detour 
electromagnetic waves around a spherical region of space, with subsequent work on arbitrarily sized and 
shaped solids by Nicolet et al. [28] . As already mentioned, the experimental validation of these theoretical 
considerations was given by Schurig et al. [16j , who used a cylindrical cloak consisting of concentric arrays 
of split ring resonators which made a copper cylinder invisible to an incident plane wave at 8.5 GHz (via 
artificial magnetic activity). However, in order to broaden the range of frequencies over which a cloaking 



^We avoid the word cloaking here as its meaning varies significantly between authors: for instance, 
Kohn et al. consider in [Q that cloaking occurs when two different media have the same boundary 
measurements, while Greenleaf et al. further require in |10) that one medium contains a cavity in 
which general objects can be hidden. The latter definition corresponds to invisibility cloaks designed by 
physicists, see e.g. [TSl US], and seems to require a singular transformation F, although whether this 
is really mathematically necessary is still unknown for n > 3. Another commonly used terminology is 
external cloaking which refers to a countable set of dipoles hidden by the presence of a resonant object, 
see [nillH]. 
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metamaterial works, one needs to explore other routes avoiding resonant elements. Huang et al. proposed 
a concentric multilayered design of cylindrical electromagnetic cloaks [53] ■ Farhat et al. further divided 
these layers in infinitely conducting sectors in the context of electromagnetics (301 • ^ very comphensive, 
yet rigorous, mathematical account of the homogenization approach to cloaking can be found in a physics 
paper by Greenleaf et al. [12], while numerical illustrations of this theory are discussed in Farhat et al. 

There are nevertheless genuinely non-singular routes towards cloaking: in 2008 Pendry and Li pro- 
posed to design an invisibility carpet so that a bump surrounded by a transformed medium mimicks the 
reflection of a flat ground, both in electromagnetic and acoustic contexts [26l|27]: Here, the geometric 
transform is a one-to-one correspondence between a segment and a curve. This proposal already led to 
theoretical |33j and experimental j34| demonstrations of three-dimensional electromagnetic carpets. An 
experimental validation of two-dimensional cloaking at visible wavelengths has even been reported for 
surface plasmon polaritons at a metal dielectric interface |35| . Also woth mentioning is the proposal of 
a non-singular cloak by Leonhardt and Tyc based on the stereographic projection of a hypersphere on a 
hypcrplane, in a way similar to what is done to design a Maxwell fisheye lens |31j . 

It is also possible to apply geometric transforms to other wave equations: However, Milton, Briane, 
and Willis |36j have shown that the elasticity equations are not invariant under coordinate transformations 
and consequently that if cloaking exists for such classes of waves, it would be of a different nature than its 
acoustic and electromagnetic counterparts. A systematic investigation of acoustic cloaking started with 
Cummer and Schurig [20j , who analyzed the two dimensional cloaking for pressure waves in a transversely 
anisotropic fluid by exploiting the analogy with transverse electromagnetic waves. Chen and Chan [23 1 124 ] 
and Cummer et al. [25] further noticed that a three-dimensional acoustic cloaking for pressure waves in 
a fluid can be envisaged since the wave equation retains its form under geometric changes. Farhat et 
al. also looked at the limitations of square acoustic cloaks for antiplane shear waves [H] . A promising 
avenue towards the realization of acoustic metamaterials was opened by Torrent and Sanchez-Dehesa 
[21] and Chen et al. [23], who independently investigated cloaking for concentric multilayers behaving as 
anisotropic fluids in the homogenization limit. Using a similar approach, Farhat et al. [22] theoretically 
demonstrated cloaking of surface liquid waves for a microstructured metallic cloak between 10 Hz and 15 
Hz, which was experimentally validated at 10 Hz. 

However, there are other physical situations whereby the mathematical model for acoustic waves is 
less tractable. For instance, Norris investigated some general types of acoustic cloaks with finite mass 
consisting of so-called pentamode materials, which display an anisotropic stiffness [37]. Also, when one 
moves to the area of coupled pressure and shear elastic waves, the isomorphism between the Navier equa- 
tions and the scalar wave equation is lost and computations become more involved. Brun, Guenneau, 
and Movchan studied a cylindrical cloak for in-plane elastic waves which is described by a rank-4 (non- 
symmetric) elasticity tensor with 2'* Cartesian entries and an isotropic density [38]. Whereas the former 
structured pentamode metamaterial |37l might already represent a technological challenge for mechanical 
engineers, the latter proposal [38] imposes even severer constraints on the material parameters. Moreover, 
the required material properties for a three-dimensional elastic cloak remain elusive thus far, as these 
would involve a rank-4 elasticity tensor with up to 3* spatially varying nonvanishing Cartesian entries. 
However, in the special case of thin-elastic plates, whose spectral properties require ad hoc numerical 
techniques [39], it has been shown that the elasticity tensor can be represented in a cylindrical basis by 
a diagonal matrix with two (spatially varying) nonvanishing entries. It seems therefore quite natural to 
start designing such a cloak before investigating the other cases. 

In this paper, we investigate the original design of singular cloaks in section 2, then turn to the analysis 
of non-singular cloaks through the blowup of a small ball in section 3 (mimetism) and the blowup of a 
line/surface in section 4 (invisibility carpets). We further show in section 5 that an anisotropic cloak can 
be approximated via a radially symmetric multilayered cloak filled with piecewise constant isotropic bulk 
modulus K and mass density p to make an object surrounded by such a coat neutral for pressure waves in 
a fluid background. Finally we consider in section 6 some cloaks with non-smooth boundaries. Section 7 
concludes the paper with some synthesis of our numerical results and some perspectives for future work. 
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2 Transformational acoustics: Spatially varying anisotropic 
densities and bulk moduli 

The key concept in developing an invisibility cloak in electromagnet ism and/or in acoustics is a simple 
one: it is the idea that sound and light see space differently. For sound/light, the concept of distance is 
modified by the acoustic/electromagnetic properties of the regions through which it travels. In geometrical 
acoustics/optics, we are used to the idea of the acoustical/optical path; when travelling an infinitesimal 
distance ds, the corresponding acoustical/optical path length is nds, where n is the local refractive 
index, equal to in the electromagnetic setting (with e the permittivity and fj, the permeability of 
the medium), and equal to \/~pJ~K (with p the fiuid density and k its compressional modulus). This 
is the first hint that it may be possible to mimic a transformation of a region of space by using an 
equivalent transformation of acoustic/electromagnetic properties. We report in Figure [l] some finite 
element computations for the scattering of a pressure wave by a rigid obstacle on its own (lower panel) 
and when it is surrounded by a singular cloak (upper panel) obtained from blowing up a point in a ball. 



2.1 Change of coordinates: Gradient behaviour. 

If one considers the geometric transform from the coordinate system {x,y,z) to {u,v^w), the relation 
between the gradient operator in both coordinate systems is given by: 



dx 

d_ 
dy 



— 



/ du dv dw \ 

dx dx dx 

du dv dw 

dy dy dy 

du dv dw 

dz dz dz 



( A \ 

du 

d_ 
dv 

d_ 



(5) 



where J^x — (J^m) ^ = d{u, v, w) /d{x, y, z) is the Jacobian of the transformation (u, v, w) 



{x,y,z). 



2.2 Transformational acoustics and weak formulation of the wave 
equation 

We now invoke the three-dimensional acoustic equation: 

2 

V. (p^^Vp) +^p==0, in f], (6) 

where $7 is a bounded domain of TZ'^ . 

Multiplying Q by a test function (j) and integrating by parts, we obtain: 



dV {yix,y,z)^-Br'^ix,y,z)p) + f dV (c^^^-^^) = (7) 



where we have assumed that the surfacic term is zero (corresponding to a Neumann integral over the 
boundary dVl). 

We now apply to equation Q the coordinate change {x, y, z) — > [u, v, w) and using Q, we deduce: 
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Figure 1 : 1 Real part of the pressure field p scattered by a rigid obstacle (Neumann boundary conditions) 
of radius 0.5m on its own (c,d) and when it is surrounded by a cloak of inner radius Ri = 0.5 m, outer 
radius R2 = \ u\ (a,b) for an incoming plane wave of wavelength A — 0.3m. (a,c) Three-dimensional 
plots, the pressure wave is propagating along the x direction (from left to right); (b,d) Two-dimensional 
plots in the (zx)-plane, the pressure wave propagates along the vertical x direction (from bottom to top). 
We note the strong backward scattering in panels (c) and (d). 



dV { (jLV(„,.,^)^.£-^JLV(„,.,^)p) det{3„,)} + J^^ dV {det{3^^)u^K'^p4>) = (8) 
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And this can be expressed in terms of J 



- I ((V,„,.„0)'' '=|3i|iv,.„,„,p) + 1 4V' {^^^'Pi) - (9) 

The new parameters p and k after the change of coordinates can be now deduced, and are given by: 



in terms of S^x O'" 



in terms of J^u 
(10) 



3 Blow-up of a small ball: Non-singular cloaking and 
anamorphism 

In 2008, Kohn and co-authors [5] have regularized the transform of Pendry et al. which blows up a point 
onto a sphere [15 . The former authors consider a small ball of radius i?o a-nd define its blow up into a 
unit ball in the following way: 



< r < iT!o, 



• Ro<r<R2, 



R2 — Ri Ri — Rq 

r + R2 

i?2 — -Ro R2 — Ro 



Importantly, this transform is continuous and piecewise smooth and this transform reduces to the identity 
on the boundary r = i?2 ■ This latter property ensures that the cloak will be impedanced matched to the 
surrounding medium, in order to avoid any reflection on the cloak outer boundary. Diatta et al. realized 
that such a non-singular approach to cloaking leads to anamorphism [421 144j . whereby a rigid object 
scatters pressure waves like another rigid object (of different size and/or shape). 
The Jacobian of the transformation is given by : 



-1 











1 











1 



with 



Ri 

Rq 

R2 — Ri 



i?2 — Ro 

Thus, the transformation matrix for the two media is given by 



{r < Ro) 



(i?o < r < R2) 













-1 








a-1 



a 



with 



Ri 
Rq 

i?2 — Ri 
i?2 — Ro 



P^R. 



Ri — Ro 
i?2 — ^0 



{r < Ro) 
{Ro<r< R2) 
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Figure 2: Real part of the pressure field p scattered by a non-singular cloak of inner 
radius 0.5m and outer radius Im for an incoming plane wave of wavelength A = 0.45m 
propagating along the vertical z-direction (from bottom to top); (a) Three-dimensional 
view for ro = 0.2m; (b) 2D view in the yz-plane for tq = 0.2m; (c) 2D view in the yz-plane 
for ro = 0.05m; (d) 2D view in the yz-plane for tq = 0.45m; We note that the larger tq, the 
smaller the scattering. In the limit case ro = 0.5m, the cloak is filled with homogeneous 
isotropic fluid. 
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We note that has two constant non-zero eigenvalues, and one spatially varying eigenvalue which 
tends to = Ri/Rq when r goes to zero. Such an invisibility cloak, which is made of two heterogeneous 
anistropic media (inside the ball r < Ri and inside the shell i?i < r < R2) is therefore non singular. 
We report in Figure [2] some finite element computations showing that when the cloak is designed from 
a ball with very small radius (c), there is no field inside the inner region of radius i?i (nearly a singular 
cloak), while for a radius Rq ~ i?i the cloak and the object inside it do not curve the wavefront of the 
incident wave (d): in that case the anisotry is very small. The case Rq = corresponds to a cloak with a 
heterogeneous anisotropic shell i?i < r < R2, and some vaccuum inside the ball r < Ri. In this way, the 
domain r < i?i is now an invisible domain for all exterior observator. However, we note that in this case, 
the first eigenvalue tends to infinity when r tends to Ri so that now has two constant eigenvalues 
and one which tends to infinity on its inner boundary r = Ri. One of two things can thus happen: one 
either considers a regular cloak with an invisibility region filled with metamaterial; or one considers a 
singular cloak with an invisibility region filled with vaccuum. These two cases have their own advantages 
and downsides: The former leaves no room for a object to hide, whereas the latter can not be achieved 
in practice. 



4 Blow-up of a line: Carpets for curved grounds 

Another type of cloak is discussed in this section: one that gives all cloaked objects the appearance 
of a flat conducting sheet in the context of electromagnetism [261 EH [33] , or a flat rigid ground in the 
context of acoustics [17]. It has the advantage that none of the parameters of the cloak is singular 
and can in fact be made isotropic using quasi-conformal mappings. It makes broadband cloaking in 
the optical frequencies possible, as chiefly demonstrated for surface plasmon polaritons at a wavelength 
of 800 nanometers (requiring a specific arrangement of about 100 dielectric pillars 200 nanometers in 
diameter) |35j . This alternative approach to invisibility cloaks takes advantage of the principle that a 
perfect conducting sheet is inherently invisible to a specifically polarized incidence when the electric field 
is perpendicular to the conducing sheet. Then if there is a bump on the ground, invisibility can be 
realized by guiding the incident wave around the bump, keeping the electric field always perpendicular 
to the curved surface, which is the function of the carpet cloak proposed by Li and Pendry. 
Let us consider the linear geometric transform : 

X = X, y z = r z + zi[x,y). (11) 

Z2 x,y 



The linear transform is expressed in a Cartesian basis as: J^x' — | ^ ^ I where a — 

dz dz 1 

Z2{x, y) i(x, y) fj-Qj-Q chain rule 

Z2Kx,y) 



1 

1 

dz dz 1_ 

dx' dy' Q 



dz z' — Z2 dzi z' — z\ dz2 dz z' — Z2 dzi z' — z\ dz2 

dx' ^'^[z2-zifdx ^\z2 ~ zif dx ' dy' ^'^ ~ zif dy ^\z2 ~ zif dy ' 

This leads to the inverse symmetric tensor T^^ which is fully described by seven non vanishing entries 
in a Cartesian basis: 

(r-i)n = (T-i),2^ , (r(-i))i3 = (T-i)3i = 

(T-1)23-(T-1)32 = -^ ,(r-l)33=(^l 

We report some finite element computations for a carpet with parabolic shape in Figure [Sj A pressure 
wave incident from above on a flat ground with a parabolic bump (lower panel) is backscattered with a 
disturbed wavefront, while the wavefront is flattened when the bump is surrounded by a carpet (upper 
panel). 
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Figure 3: Lower panel: Pressure wave of wavelength A — O.Sm incident from above on a flat (rigid) 
ground with a bump described by the surface z\(x, y) = \J \ — — — 1.802776; Upper panel: Pressure 
wave of wavelength A = 0.3m incident from above on a carpet of parabolic shape with inner and outer 
boundaries Z\(x^y)) and Z2{x,y) = y/l — x"^ — y'^ — 1/2 surrounding the same bump as aforementioned; 
(a) and (b) are three-dimensional plots, while (c) and (d) are two-dimensional plots along the vertical 
slices. 
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5 Approximating ideal cloaks: Structured cloaks 



Transformational acoustics is a nice way to imagine new kinds of metafluids, but it is however important 
to be able to fabricate these with materials/fluids at hand. To achieve this practical goal, one possible 
route is to use homogenization techniques which have been developed over the past 40 years by applied 
mathematicians [40l Hi] . Such multiscale techniques make effective medium theories rigorous and further 
allow for more flexibility in the geometrical and material parameters: one can homogenized materials 
with small inclusions of any shapes in two or three space dimensions, and these can be anisotropic. 



5.1 Homogenization of the acoustic equation 

In our case, the pressure field is solution of the following governing equation: 



V- (A„Vp) + ^C.,P = 0, (13) 



c 



inside the heterogeneous isotropic cloak ilf, where 

C^ = k"^(-) and A^=/9"^(-). 

Furthermore, = kq/poj where kq is the bulk modulus of the fluid and po its density. 

When the acoustic wave penetrates the structured cloak it undergoes fast periodic oscillations. To 
filter these oscillations, we consider an asymptotic expansion of the associated pressure p solution of ([T3| 
in terms of a macroscopic (or slow) variable x — {r,9) and a microscopic (or fast) variable x^ ~ (^'^)' 
where 77 is a small positive real parameter. 

The homogenization of this elliptic equation can be derived by considering the following Ansatz 



Vxe%, p,(x)=^77>«(x,x,), (14) 
=0 



where Vx (z flf, p'-^H^; ') is 1-periodic along r. Note that we evenly divide fi/ (i?i < r < R2, < 6 < 2tt) 
into a large number of thin curved layers of radial length {R2 — Ri)/ri, but the spectral parameter /3q in 
Eq. ( 13 1 remains fixed (so is the wave frequen cy) . 

Rescaling the differential operator in Eq. 



1 d 

13 1 accordingly as V = Vx + - — , and collecting the 



terms sitting in front of the same powers of 77 (see e.g. [101 HI]), we obtain 



77 dr 



V • (AVp) + — < c > P - , (15) 



which is the homogenized acoustic equation where < C >— Jq ^{r)dr and with A a homogenized 
rank-2 diagonal tensor which has the physical dimensions of a homogenized anisotropic density p — 

Dia.g{pr, pe,P4,) given by 

A = Diag(< < \>,< X>) = p-^ . (16) 

We note that if the cloak consists of an alternation of two homogeneous isotropic layers of thicknesses 
d-A and ds and bulk moduli ka, kb and densities pA and pB, '^g have 

1 1 f 1 v \ PA + ripB KA + VKb 

— = 7-^ — \ , Pe^ P^ = — — , K — — , 

Pr l + r]\pA Pb J 1 + f? 1 + ?/ 

where rj = ds/dA is the ratio of thicknesses for layers A and B and + d_B — 1- 

We now note that the coordinate transformation r' = Ri + r ^^^^^ can compress the region r < R2 
into the shell Ri < r < R2, provided that the cloak is described by the following density and bulk 
modulus 

R2-Ri( r V (R2-Ri\ fR2~RiYf r Y 
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where Ri and i?2 are the interior and the exterior radii of the elastic coat of thickness h. 

To mimic these fluid parameters, we proceed in two steps, following |23| : we first approximate the ideal 
cloak by a multi-layered cloak with M anisotropic homogeneous concentric layers. We then approximate 
each layer i, i — 1, .., M by N thin isotropic layers through the homogenization process described above. 
This means the overall number NM of isotropic layers can be fairly large. 



5.2 Numerical results for a multi-layered cloaking device 

To check the functionality of the multi-layered acoustic cloak, we have performed some finite element 
simulations. Maps of the pressure wave at wavelength 0.25m propagating at a time instant are repre- 
sented by the real part of the complex amplitude p from the left in Fig. |4] and is scattered by a rigid 
circular obstacle (left panel a) and by the layered cloak (right panel b). When it is surrounded by the 
heterogeneous cloak, which consists of a circular coating of radii Ri = 0.2 and R2 = 0.4 meters, both 
forward and backward scattering nearly vanish, with interestingly a phase shift between the acoustic 
wave propagating in homogeneous fluid and the wave bent by the cloak (the acoustic path is obviously 
much different). We note that this phase shift cannot be avoided in electromagnetic metamaterial cloaks, 
as otherwise the wave propagating within the cloak would have to travel faster than the speed of light 
in vacuum. The acoustic parameters of the proposed layered cloak are characterized by a scalar bulk 
modulus K and a spatially varying rank 2 tensor p given by (17), see Figure H] for the plots of the pressure 



field scattered by such a cloak and Figure [5] for the parameters of the cloak. Moreover, we numerically 
checked that the cloak is broadband in nature. 

To further demonstrate the validity of our approach to broadband cloaking, wc now propose an 
original design for the non-singular cloak proposed by Kohn et al. [61 (i.e. deduced from the blow up 
of a small ball instead of a point) with 20 concentric layers of homogeneous fluids, see Fig. |6] We note 
that the constraint on the fluid parameters is relaxed (the lower bound for the density is higher). We 
hope this simplified set of material parameters will pave the way towards an experimental realization of 
a near-ideal acoustic cloak in the near future. 
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6 Mimesis with Polyhedron and Star 



One of the specific features of polyhedra and stars is their symmetries, allowing to split the cloak into 
sub-regions with (smooth) linear boundaries. Such sub-regions are diffeomorphic via a linear action of 
elements of a symmetry group. This allows to simplify the design of such cloaks as we study one fixed sub- 
domain, and deduce the properties of the other remaining ones using the diffeomorphisms mapping it to 
the other ones, which are parts of a symmetry group acting on the whole domain. For the two-dimensional 
case, we refer the reader to |33]. 



6.1 Design of three-dimensional cloaks with corners 

For this construction, we use spherical coordinates (p, 9, 0) related to cartesian coordinates by 

p — ^/x^ + + z^, 

' = arctan(y/a;), (18) 
> = arccos( , ^ ) 

where p, 9 and (j) belong to the appropriate subsets of M in accordance to the geometry of the cloak. Or 
equivalently, we have 

X ~ pcos(9) sin(0), 

y = psm{9) sin((/)), (19) 
z — pcos{(t>) 

We denote by 3xp the Jacobian of the latter change of coordinates. That is, 

(cos{9) s'm{(f)) — sm{9) sin((/)) cos{9) cos{(j))\ 
sm{9) sin(0) cos{9) sm{(f>) sm{9) cos(0) (20) 
cos(0) — sin(0) J 

To express the transformation (p, 9, <j)) i-> [p'{p^ 9, (j)), 9'{p, 9, (j)), (j>'{p, 9, we consider three noncollinear 
points (a;i,i,?/i,i, Zj_i), {xi^2, yi,2, Zi,2), (a;i,3, 2/i,3, z^s), on each Si of the boundary surfaces 5*0, ^i, 5*2 of the 
sub-domain, then we use 

^'-^i + tS (P-^O), 
9' = 9, 

where Ri are now smooth functions R{9, (p) of 9 and (j)) 

j;^ ^ 12 f 2 1 1 

* ai cos{9) sin{4>) + hi sin(0) sin((/)) + Ci cos{4>) ' 

and ai, hi, Ci, di are constants 

a-i = {yt,2 - yi.i)izi^3 - Zi.i) - {zt,2 - Zi^i){yi^3 - yi^i), 

bi = iZi^2 - Zi^i){xi^3 - Xis) - iXi^2 - Xi^l){zi^3 - Z^^l), a = {Xi^2 - Xi^i){yi^3 - yis) - iyi.2 - yi,l)iXi^3 - Xi^i), 

di = {yi,2Zi.3 — yi,2Zi^l — ViAZi^S — 2^1,22/1,3 + ^i,22/i,l + 3)2^1,1 + (^i, 22^1, 3 ~ Zi 2Xi^i — 2i, 1X^,3 — Xi,2^i,3 + 2;^, 22^,1 

+ 2;i,iZi,3)j/i,i -I- (xi,22/i,3 - a;i,2?/i,i - x^^iyi^s - yi,2Xi,3 + yi,2Xr,i + 2/i, 13:^,3)2:^,1. (22) 



Now if we set a 



R2-B1 
R2-B.0 ' 



_ dRi _ dijai sin{9) sin(0) - hj cos(6>) sin(0)) 
^' 09 {tti cos{9) sin((/)) -I- hi sm{9) sm{(j)) + a cos{(j)))^ ' 



Ri 



dRi —diCi sin((^) 



d4> (flj cos(6') sin((/)) -I- hi sin(6') sin((/i) + c; cos{4>))^ ' 
da 
09' 
Op' 



Oa Oa I Op' ^ / n \ n 

1)9' ^'^ ''^ ~09 ^ ^ ^ - ^0) - oiRo^i 



Ri,4> + oi^{p — Rq) ~ ctRo,(t,, 
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The Jacobian 3p>p of the above transformation {p, 9, (f) i-> {p'{p, 9, (f)),9'{p, 9, (j)), (j)'{p, 9, (f>)), is 



Jp'p — 



(24) 



However, for implementation, one may need to get back to Cartesian coordinate system. The following 
diagram gives the picture as to how to compute the resulting Jacobian J^'x of change of coordinates, in 
Cartesian coordinates. 



{x,y,z) 



{x',y',z') 



' p'x' 



So it now reads 3x'x — J x' p' p' p-f px 

^ 3^diag{l,p',p')3p'pdiag{l,^,pJ^ . We are now ready to 
deduce the material properties by the usual formula T^^ — Jx'x^x'x/ det(Ja;'a;). Note that 

Jx'x^l'x = J4,diag{l, p' , p')3p>pdiag{l, ^, ^)3^'^3^'^ diag{l, ^, ^)3'^,pdiag{l, p' , p')3l (25) 

and using JTJ^ = diag{l, sin^(^), 1) it becomes 



^x'x^l'x = 3^diag{l,p',p')3p,pdiag{l, ^, ^)diag{l, . ^ , l)diag{l, -, -)3Z diag{l, p' , p')3l (26) 

p p sm (0) ^ ^ 



1 



1 1, 

P' P' 



so that 3x'x3'^,^ = 3^diag{l,p',p')3p,pdiag{l, p^;^^, j2)3j,pdiag{l, p' , p')3l. 
This way, the material coefficients will read in Cartesian coordinates 



det{3x'x) 



Pe 

p^ sin-^{0) 

Pe 
sin^ 



P<t. 



Pe 



p2 sin^(0) p' 

1 

j,) 



(27) 



6.2 Numerical illustration for a Regular Polyhedron (isocahe- 
dron) and its stellated counterpart (6-point star) 

We first consider a linear geometric transformation that maps the point of coordinates (0,0,0) on an 
isocahedron Ii of edges of length ei — 0.2m and that maps Ii on an isocahedron I2 of edges of lengths 
62 = 0.4m, (ei < 62). One can see in panels (a) and (b) of Figure [T] that the pressure field radiated 
by an acoustic source is smoothly bent around the inner region of the cloak (the amplitude of the field 
indeed vanishes in this region) . When the number of edges increases, the cloak's boundary approximates 
to certain extent that of a spherical cloak, and the scattering is reduced, see [13] for a detailed analysis 
in the two-dimensional case. We then consider a linear geometric transformation that maps the points 
of coordinates (0, 0, 0) on a 6-point star 5*1 of edges of length ei = 0.12m and that maps Ii on a 6-point 
star 5*2 of edge of lengths 62 = 0.45m, (ei < 62). One can see in panels (c) and (d) of Figure [t] that the 
pressure field radiated by an acoustic source is smoothly bent around the inner region of the star-shaped 
cloak. When one increases the number of edges (hence of corners, which become sharper and sharper), 
the cloak's boundary becomes more singular, and the scattering worsens, as shown using an asymptotic 
approach in the two-dimensional case in [43] . Such a behaviour is antagonistic with that of the polyhedra 
cloaks. 
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7 Concluding remarks 



In this paper, we first reviewed the chronology of mathematical results in inverse problems which pre- 
date the 2006 proposals of Pendry et al. [B] and Leonhardt [9j for an invisibility cloak. We trace 
these results back to 1984, when Kohn and Vogelius noted the arbitrariness of the Dirichlet-to-Neumann 
map for conductivity problems involving matrix valued coefficients [1], with subsequent works by other 
mathematicians including Lassas, Uhlmann and Greenleaf who published the first cloaking paper in 2003 
|32j . We reviewed the formulae underlying transformation based acoustic cloaks and performed some 
finite element numerical simulations for singular and non-singular cloaking effects in three-dimensional 
settings. It has been recently shown that acoustic cloaking shells are possible by means of multi-layered 
structures consistiing of two types of acoustic isotropic materials, whose acoustic parameters should be 
radially varying |21| . Here, we proposed a design of a broadband multi-layered cloak with a large number 
of thin homogeneous isotropic layers, based on a geometric transform for a non-singular cloak proposed 
by Kohn et al. _6J. For this, we first rederived the homogenized acoustic 3D equation using a multiscale 
asymptotic approach. We found that the homogenized acoustic parameters are described by a rank-2 
tensor (a generalized density) and a scalar bulk modulus both of which are functions of the radius. We 
then performed numerical computations based on the finite element method which proved that a rigid 
obstacle surrounded by a coating consisting of 20 concentric layers alternating two types of isotropic bulk 
moduli and mass densities is neutral (vanishing backward and forward scattering) for the acoustic plane 
waves, over a finite interval of frequencies (in Hertz). The cloak can also mimick the scattering of a 
small rigid obstacle. We finally investigated 3D acoustic cloaks with an irregular star-shape, in order to 
exemplify the versatility of transformation acoustics. We emphasize that our numerical results can have 
interesting applications in sonar camouflaging in the context of anamorphism (mirage) effects. 
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Figure 4: Real part of the pressure field p scattered by a rigid obstacle of radius Ri = 0.2m surrounded 
by a non-singular multi-layered cloak of inner radius Ri and outer radius i?2 — 0.4m for an incoming plane 
wave of wavelength A = 0.25m. These cloaks scatter waves like a rigid obstacle of radius ro (mimetism). 
Upper panel: Cloak scattering like a rigid sphere of radius — 0.05m , the ranges of parameters are 
Pi/Po = [0.0243,3.4636] and k^/ko — 0.5714; Lower panel: Cloak scattering like a rigid sphere of radius 
To = 0.15m, with ranges of parameters Pi/po — [0.2601,2.2229]; Ki/Ko = 0.8. We note that the density 
takes less extreme values for tq — 0.15m. (a,c) Three-dimensional plot (the pressure wave is incident 
from the left); (b,d) Two dimensional plot in the vertical zx-plane (the pressure wave is incident from 
below) . 
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Figure 5: Variation profil of the reduced parameter Prj Pq versus the thikness of the cloak (red curve) 
and a layered approach (blue staircase curve). Here, we consider the blowup of a small ball of radius 
ro = 0.05m (a) and a small ball of radius 0.15m (b). Importantly, the lower bound for the density in (b) 
is about five times larger, while the upper bound is similar to that in (a). It should be also noted that 
the graph in (b) is the restriction of the graph in Fig. |6]to the interval r e [0.2m, 0.4m]. 
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Figure 6: Real part of the pressure field p scattered by a non-singular multi-layered cloak of inner radius 
Ri = 0.2in and outer radius i?2 = 0.4ni for an incoming plane wave of wavelength A = 0.25m. We consider 
the blowup of a small ball of radius rp = 0.15m. (a) Three-dimensional plot (the pressure wave is incident 
from the left); (b) Two dimensional plot in the zx-plane (the pressure wave is incident from below); (c) 
The variation of the bulk modulus versus the radiusq is given by the dashed black curve, while the 
variation of the density (solid red curve) is approximated by the blue stair-case curve. Importantly, the 
cloak is filled with a homogeneous isotropic fluid of density Pr/ Po = 0.75 and bulk modulus k^/kq = 1.333. 
Importantly, the density and bulk modulus display a jump at the inner boundary of the cloak i?i and they 
are both constant in the core i.e. for r < Ri (isotropic homogeneous fluid). It is therefore counter- intuitive 
that the system core and cloak is invisible (while the core and cloak taken separately arc not). 
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Figure 7: Upper panel: Scattering by an isocahedron cloak in presence of a point source; (a) Isocahe- 
dron cloak 3d-view; (b) Slice in the (x, z)-plane for isocahedron with edges ei = 0.2m and 62 = 0.4m at 
the frequencie A = 0.3m. Lower panel: (c) Geometry of the 6-point star cloak (3d-view); (d) Geometry of 
the (x,y)-plane for a star with edges ei — 0.12m and 62 = 0.45m at the frequency A — 0.3m; Lower-Left: 
Scattering for a point source in presence of the star cloak (3d- view); Lower- Right: Scattering for a point 
source in presence of the star cloak (2d- view). 
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